Phase Separation and an upper bound for A for Fermi fluids in the unitary regime 
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An upper bound is derived for A for a cold dilute fluid of equal amounts of two species of fermion 
in the unitary regime kja — > oo (where kf is the Fermi momentum and a the scattering length, 
and A is a pairing energy: the difference in energy per particle between adding to the system a 
macroscopic number (but infinitesimal fraction) of particles of one species compared to adding equal 
numbers of both. The bound is S < § (2(2£) 2/5 - (2£)) where £ = e/e FG , S = 2A/e FG ; e is the 
energy per particle and efg is the energy per particle of a noninteracting Fermi gas. If the bound 
is saturated, then systems with unequal densities of the two species will separate spatially into a 
superfluid phase with equal numbers of the two species and a normal phase with the excess. If 
the bound is not saturated then A is the usual superfluid gap. If the superfluid gap exceeds the 
maximum allowed by the inequality phase separation occurs. 



43 

o 



I 

C 

O 
o 



> 
o 

oo 
O 

o 
in 
o 

c3 



-a 
c 

o 
o 



X 



PACS numbers: 

During the past several years there has been consider- 
able theoretical interest in studies of cold, dilute, Fermi 
systems with equal densities of two strongly coupled 
species (HUE S 11 S H H E3 In this context species 
refers both to spin states and to internal quantum num- 
bers. There is an idealized version of this problem in 
which the following conditions are met: the particles of 
each species are of equal mass and 1/a <C n 1 l z < 1/ro 
and T <C n 1 / 3 where tq is the typical distance scale of 
the interaction a is the scattering length for free parti- 
cle scattering between the two species, n is the density 
of each species and T is the temperature and units have 
been chosen with % = 1. In the extreme limit of this 



situation where 



oo, Aa 



oo, n 1 / 3 /A 







and Tu -1 / 3 — > 0, (where A is the typical momentum 
scale characterizing the interaction) there is only a sin- 
gle momentum scale in the problem, namely n 1 / 3 . This 
is sometimes called the unitary regime. It is conve- 
nient to re-express this in terms of a nominal Fermi mo- 
mentum kf = (67r 2 n) 1 / 3 . Thus all physical observables 
in the problem can be expressed as appropriate pow- 
ers of this scale times appropriate constants. For ex- 
ample the average energy per particle can be written as 

e = £ 1 = £ £ fg where £ is a universal constant (epc 
is the energy density on a noninteracting Fermi gas). A 
pairing energy parameter which gives the difference in 
energy per particle between adding a macroscopic num- 
ber (but infinitesimal fraction) of particles of one species 
as compared to adding equal numbers of both to the sys- 
tem with equal particle numbers can similarly by given 

as 2 A = <5|tt^ = SepG- The connection between A and 

5 2m r ^ 

the usual superfluid gap is somewhat subtle and will be 
discussed below. 



This problem is of interest in part due to the univer- 
sality of the behavior. The coefficients such as £ and 
6 apply to all problems in this regime regardless of the 
microscopic details of the problem. The problem also 
is of theoretical interest in that it represents the exact 



intermediate limit between two weakly coupled regimes 
with kfci <C 1. The true weak-coupling regime between 
fcrmions has kfa<0 and is a BCS superfluid; the regime 
with kfa small and positive corresponds to weakly cou- 
pled molecules in a BEC 11]. The problem is also chal- 
lenging there appears to be no simple analytical method 
to compute the universal coefficients. 

The problem is relevant to physical systems of inter- 
est. In nuclear physics the problem of low density neu- 
tron matter can be caricatured by such a system: the 
two species are the two spin states of the neutron; the 
s-wave scattering length between spin up and spin down 
neutrons is much larger than the characteristic range of 
the nucleon-nucleon force 0, Hi! flfij ] . The problem has 
become of importance at the interface between atomic 
and condensed matter since the scattering length be- 
tween atoms in particular m states can be tuned via al- 
tering an external magnetic held. The scattering length 
diverges at a Feshbach resonance. There has been in- 
tense experimental work on pairing in fluids of trapped 
Fermionic atoms and the transition from the BCS to BEC 
regime 12]. Of course, the trap itself can play an impor- 
tant dynamical role in the problem and significant theo- 
retical effort has gone into describing the role of the trap 
which adds a spatial dependence to the problem [ljj. The 
present paper will focus be on the ideal case where the 
particles are visualized as being contained in a large box. 

The fact that no direct analytical computations of 
the relevant dimensionless parameters exists means other 
methods must be found to learn something about these 
dimensionless parameters. One strategy is attempt to 
extract them numerically 0, 0, Q, @- A possible diffi- 
culty with such an approach is that a priori estimates 
of the errors may be difficult to obtain in a reliable way. 
Thus, a constraint based on reliable analytical methods 
is potentially quite useful. One possible idea is to see 
whether the coefficients £ and S can be related to each 
other analytically. This is is also a formidable challenge 
for which no rigorous answer is known. However, as will 



2 



be discussed in this paper it is possible to give a rigorous 
upper bound on S for any assumed value for £: 



5 < 



(2(2£) 2/5 - (20) 



(1) 



As will be discussed below, A need not be the superfluid 
gap and the superfluid gap can exceed this bound; if it 
does one predicts an interesting phenomenon in the case 
were the densities of the two species are unequal. 

To derive this bound, consider a generalization of the 
problem to the case where the two species (denoted a 
and b) have different number densities: the only rele- 

1/3 

vant quantities with dimension of inverse lengths n a 
are n^ 3 . One completely general way to parameterize 
the ground state energy density of this system subject to 
the constraint of fixed density of the two species consis- 
tent with the correct dimensional scaling is 



£(n a ,n b ) = an 5 J 3 f(n b /n a ) 



(2) 



where / is a universal function which only depends on 
the ratio of the number densities and a is a constant with 



dimension of mass -1 . At the point rib = the system is 
a noninteracting Fermi gas of species a (by hypothesis 
the only relevant interactions are for species a and b to 
interact with each other). With out loss of generality 
one can fix /(0) to be unity and this in turn fixes a to 
its Fermi gas value. The parametrization in Eq. J2J) is 
very natural if one envisions starting with a Fermi gas 
of species a and slowly adding in particles of species b. 
Comparing Eq. (J2J) with the definition of £, one sees that 
£ = /(l)/2. The factor of 1/2 in this relation reflects 
the fact that at x = 1 the two species contribute equally; 
were they non-interacting /(l) = 2 and £ is defined as 
the fraction relative to the noninteracting case. 

The key thermodynamic consideration to derive a 
bound is the possibility of phase separation. This con- 
strains the energy density as a function of the densities. 
In particular, if we consider the energy density £(n a ,rib) 
at two different pairs of number densities for the two 
species, (na,n^) and (n« ,n b ), then the average of 
the energy densities at these two number densities cannot 
exceed the energy density at the average number density: 



J 



£ (ni 1} , n( l) ) + £ (ni 2) , n { b 2) ) ^ „ ( + n£ 



> £ 



(2) n^+nf' 



(3) 



r 



if Eq. were false, it would be possible for a sys- 
tem with a fixed but large volume and number densi- 

4.- fn^+n™ '4 1) +«< 2) \ , , ., , ,. ... 

ties I — — 2 , — - — 2 I to lower its energy by dividing 

the volume into two equal regions with different phases: 
one with number densities (n^p, n^) and the other with 

(«i 2 U 2) ). 

Equation 13 implies that in regions where £(n a ,rib) is 
continuous its curvature in any direction in the n a ,rib 
plane is positive: 

d 2 £ 

a — A? ^ ( 4 ) 

OJli , Oflj 



for all unit vectors n where i,j can assume the value of 
a or b and summation of over i and j is implicit. 

In order for Eq. Q to hold for any unit vector n, the 
matrix 



K(n a , rib) 



d 2 £ 

'In " 



d 2 £ 
dn a dnt 



d 2 £ 
dn'r 



(5) 



must have only nonnegative eigenvalues; thus det(K) > 
0. Inserting the parametrization of Eq. @ into the def- 
inition of K and imposing a positive determinant yields 



a constraint on the curvature of the function / 

2f\x) 



f"(x) > 



5/(s) 



(0) 



We know /(0) = I and /(I) = 2£ and that / contin- 
uously connects these with its curvature constrained by 
Eq. 10. Consider the curve which obeys these boundary 
conditions and saturates inequality © at all points in 
between. Define that curve as / max : 



/max(^) 

with 



2/max ( x ) 



where = 2£ 



5/*max(*^) 
/max(0) = l /max(l) 



(7) 



The differential equation in Eq. Q can easily be solved 
subject to the boundary conditions. There is only one 
real solution: 



/max - (l + ((2£) 3/5 - l) x) 



5/3 



(8) 



The differential equation for / max was derived by con- 
sidering a path associated with varying the densities n a 
and rib which always is in the direction where the sec- 
ond derivative of £ is zero: thus the derivatives of £ with 
respect to n a (or n b ) (i.e. the chemical potentials) are 



3 



constants along the path . Therefore / max represents a 
situation in which the system for x = and x = 1 and at 
all points in between are at the same chemical potential 
(but different total density). This is precisely the condi- 
tion for phase separation: at < x < 1 a fraction r of 
the particles are in the superfluid phase with a density 
of n s and a fraction 1 — r are in the normal phase with 
a density n n . It is a simple exercise of matching chemi- 
cal potentials and densities to show that in such a phase 
separated system: 



i a (l + x((20 3/5 -l)) 



(2£) 3/5 

x(2£_) 3/5 
1 + x((20 3 / 5 - 1) 



(9) 



where n a is average density of type a over the entire sys- 
tem. Note that eqs. (jHJ and © can also be derived by 
assuming at the outset two phases and then varying their 
densities and fractions subject to the constraints of fixed 
average density and fixed ratio of total number of the two 
species. 

The phase separated configuration with fixed x and 



n a has a known energy. The actual minimum energy 
configuration is either this energy or below it so / ma x 
serves as an upper bound for /: 



f(x) < /max (x) 



(10) 



a homogeneous phase violating this condition is ener- 
getically unstable against phase separation and thus the 
ground state is phase separated; the upper bound is sat- 
urated if phase separation occurs. 

It is worth observing that the the preceding analysis 
is valid only for rib/n a < 1. However, the regime x > 1 
can easily be studied as it corresponds to more of species 
b than a. For rib > n a one can use the previous analysis 
with b and a switched: 



£(n a ,n b ) = anl /3 f(n a /n b ) 



(11) 



with / the same function as above. 

One can determine A from f(x). For an ordinary su- 
perfluid A is the gap. The gap represents the amount 
of energy saved by pairing: A is the difference in energy 
per particle gained by adding particles of one species type 
(say type a) to a system of equal particle number as com- 
pared to the energy of adding equal numbers of a or b: 



2A 



(E(N + 2M, N) - E(N, N)) - (E(N + M,N + M) - E(N, N)) 

M 



(12) 



r 



where E(N a , Nb) is the total energy and N a ( Nb) is the 
total number of particles of species a (b) and M <C N is 
the number of particles of each type added. Going to the 
thermodynamic limit gives A as the discontinuity of the 
derivative of £ with respect to the density of one of the 
species: 



2A 




(nb—na-\-e) 



d£ 

dn h 



(13) 



Using the general parameterizations of £ of Eqs. @ and 
(fTT|> yields: 



2A 



;.2 
2m 



€ 2 



6/'(l) 
5£ 



6/'(l) 



(14) 

where the second form follows since 6 = A/ ( ff^Y 

Equation l|14fl provides the basis for inequality Q). 
Note that / < / ma x in the interval from zero to unity 
and that by construction f(l) = f max (l). This is possible 
only if f'(l) > / max (l). Thus, 8 < £ (f - 2/ max (l)/C)) ■ 
Using the explicit form of / max from Eq. © immediately 
yields inequality 



The interpretation of inequality Q is subtle. A is the 
usual superfluid gap A$f in the case where the system 
does not phase separate for unequal numbers (i.e. the 
inequality is not saturated). If there is only one possible 
phase when particles are added, then 2A must simply 
represent the pairing energy for this phase. In the case 
where inequality is saturated, however, this is not the 
case. Although A retains the definition given above, it 
should not be interpreted as A$f ; if a mixed phase is 
energetically preferred, A represents the amount of en- 
ergy per particle to add particles of one species into a 
normal phase which forms in equilibrium with the super- 
fluid phase jrj|- The distinction is the following: A$f 
is the energy per particle to a single particle to a sys- 
tem with equal numbers of the two species; A represents 
the energy cost per particle when adding an large num- 
ber (but infinitesimal fraction) of particles of one species. 
Clearly, Agi? > A ; either phase separation happens and 
the two are equal or it does not and it is energetically 
cheaper to add unpaired particles in a new phase . 

To summarize: inequality |£Q)always holds with A de- 
fined as above. If in addition the system is known not 
phase separate at unequal particle numbers then i) the 
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inequality is not such saturated and ii) A = Asf- 

An important corollary of this analysis is that if Ssf > 
f (2(2£) 2/5 - (2f)) (where S SF is the analog of 5 for the 
superfluid gap) phase separation must occur for x < 1 
and if S SF < § (2(2£) 2 / 5 - (2£)) then the type of phase 
separation considered here (into a fully paired and fully 
unpaired phases) does not occur. 

Given one nontrivial physical assumption, it is possi- 
ble to make a much stronger connection between £ and 
S than inequality JQ). The dynamical assumption is that 
for n a 7^ rib, the system does separate spatially into two 
phases: a superfluid phase (which prefers to have equal 
numbers of the two species) and a normal phase with the 
remainder. This possibility was explored in an intrig uing 
recent paper by Bedaque, Caldas and Rupak (BCR)[l4j. 
The BCR paper argued on the basis of a generalized BCS 
ansatz that that such phase separation occurs. The anal- 
ysis of BCR was aimed at a broader class of problems 
than the strongly coupled limit of kta ^> 1; indeed, the 
detailed analysis is only strictly legitimate in the case of 
small BCS gaps and hence weak coupling. Thus, it is an 
open question as to whether such phase separation occurs 
at strong coupling and asymmetric systems. However, 
BCR argued that it was plausible that their conclusion 
holds even away from the weak coupling limit. 

Note, that if there is phase separation as suggested 
by BCR, the inequalities in and 1|10|) must be sat- 
urated. Since the chemical potentials for both species 
are constant at all points along the path a region with 
rib = is in chemical equilibrium a region with rib = n a ; 
this is necessary and sufficient for phase separation of the 
BCR type. Thus, the BCR assumption implies that for 
> x > 1, f{x) = (1 + ((2£) 3/5 - 1) x) 5/3 and thus 

<5=^(2(2£) 2 / 5 -(20) . (15) 

Clearly it is important to establish whether phase sep- 
aration occurs for x ^ 1. As discussed above, this can be 
immediately answered if one knows £ and Ssf (the analog 
of S for the superfluid gap). Estimates of £ and Ssf have 
been obtained numerically usingMonte Carlo methods 
for finite but large systems 0, 13 ■ The most recent 
extracted value [7j of £ is approximately .42 ± .01 which 
implies that the dividing line between whether phase sep- 
aration occurs or not is approximately Ssf = 1.70 with 
fairly small numerical uncertainty. The extracted value 
of for Ssf is 1-68 ± .1. Unfortunately this is not accurate 
enough to determine whether phase separation occurs. 
The numerical simulations in ref. Q for x < 1, are ener- 
getically consistent with phase separation. The present 
analysis implies that if these numerical simulations re- 
liable, then either S SF does exceed f (2(2£) 2 / 5 - (2f)) 
(presumably by a small amount) or that phase separation 
does not occur for x < 1 with the energy just slightly 
below the phase separated energy. Thus, the analysis 



here provides a highly nontrivial constraint on the nu- 
merics. Finally, it is worth noting that it is surprising 
just how close Ssf is to the critical value for phase sepa- 
ration: present numerical simulations do not rule out the 
intriguing possibility that they are exactly equal. 
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